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ABSTRACT 

In this paper, we demonstrate the equivalence of a scalar input system 
X =J3^ + bu, for which the eigenvalues of the generator coincide with the 
roots of the entire function 


^ (o(T-0,) 

> V (oT 1 

p(to) = e + Sj^e + 


+ a + 
m 


fa(0)e“^'’^"®^de, 

0 


with the controlled scalar functional equation 


T 

y(t) + a,y(t-0J + ••• + a y(t-T) + / a(0)y(t-0)d6 = u(t), 

1 1 m 0 

The theory of nonharmonic Fourier series is then employed to Investigate the 
placement of eigenvalues in the closed-loop system with continuous state 
feedback. 
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1. Introduction 

Let be a generator of a strongly continuous group of operators 
{^t) : t e ^ on a Hilbert space and suppose that the spectrum of 
consists of an infinite sequence of simple eigenvalues which forms the 

zero set of an entire function having the form 


p((o) = e + a^e 


0)T . 4. /ON 

+ ••• a + j a(9)e d0* 

“ 0 


( 1 ) 


In this paper, we will analyze in detail the transformation which carries the 
scalar input system 


x(t) =-Q^x(t) + bu(t), 


( 2 ) 


to the scalar functional equation 


T 

y(t) + a.y(t-6,) + ••• + a y(t-T) + / a(0)y(t-0)d9 = u(t). (3) 

X i m ^ 


The latter constitutes the causal control canonical form for the pair 
The adjective "causal” is used here because a feedback law of the form 


u(t) = (i5,x(t)) 


(4) 


will be shown to transform to a feedback law of the form 

T 

u(t) = / g(e)y(t-6)d0, (5) 

0 


i*e,, the input u depends only upon past values of y* The expression 
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“control canonical form" refers to the fact that the above-mentioned 
transformation ;has a structure quite similar to its well-known finite 
dimensional counterpart as well as to the fact that the effect of linear 
feedback (4) in the original system (2) can be readily analyzed in terms of 
the system (3) with corresponding feedback (5)* 

Before proceeding, we present an example of a system which possesses the 
structure indicated above. Consider the linear hyperbolic system in 
characteristic normal form [3]: 



with boundary conditions 

w“(0,t) = DQw'*‘(0,t) ; w'^(l,t) =D^w“(l,t). 

Here, 

A (x) = dlag(Xj(x),***,X^(x)], 

A’^(x) = dlag(X^^j (x) , • • • ,X 2 ^(x) ) , 

are diagonal n^n matrices with 

X (x) < • * • < X (x) < 0 < X (x) < • • • < X (x) . 
1 n n+l zn 

Also, 

A (x) Aj^ 2 (x) 

A(x) = 

A2j(x) A^(x) 
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and Dq and Dj are both nxn matrices which determine the manner in which 
"information" is reflected at the boundaries x = 0 and x = 1. Such systems 
arise in the study of counterflow heat exchangers [9] and in the study of 
multiconductor transmission lines [6]. 

If Aj 2 = ^21 ^ off-diagonal components of both A“ and A"^ 

vanish identically, then the characteristic frequencies are the roots of 

irCm) = det [exp(i)JA_j_ + M^) - Djexp(mA_ + M_)Dq], 


where 

1 + -1 

A = / (A*(x) ^dx, 

* 0 

and 

1 -1 4 . 

M = / (A (x) ^A*(x)dx. 

* 0 * 

The roots of tt will coincide with those of a function p having the form 
displayed in (1): 


ir(a)) 


wa/' o/T 
e [e 


aie 


o)(T-0j) 


+ 


+ a 

m 


) 


= e • p(w). 

In the case where A has non-zero off-diagonal components, one can show [19], 
under the assumption that both Dq and Dj^ are invertible, that p remains 
the same except for the addition of an integral term as in (1) with 
a(-) £ L^(0,T). 

In the above example, one may consider control problems wherein one 
introduces a forcing term proportional to a scalar function u(t) in the 
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differential equations (distributed control): 


9w 

It 


A(x) + A(x)w + b(x)u(t), 


( 6 ) 


or in the boundary conditions (point control): e.g., 


w'*'(l,T) = Djw"(l,T) + bu(t). 


(7) 


Either way, one is led to consider a system of the form 


x(t) + bu(t), 

where the element 6, in the former case, lies in the state space but, in 
the latter case, must be interpreted as a distribution* We will discuss this 
point in greater detail in Section 2* 

Of central Interest in this paper is ^the spectral synthesis problem: 
given a set of "desired** eigenvalues construct a feedback law 

of the form (4) such that the eigenvalues of the closed-loop system 

X + b( ^,x) = C^+ b <8> )x, 

coincide precisely with Russell [14] has carried out a study of this 

question in the case corresponding to ra = 1 above for a class of linear 
hyperbolic systems consisting of a pair of equations. In his study, the 
theory of nonharmonic Fourier series is used to study the canonical equation 
(3) and conclude that any sequence which 
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{(\ - \)/\l ^ (8) 

where {bj^} is the sequence of expansion coefficients of the control 
distribution element b with respect to the eigenfunctions of _q/, may be 
synthesized by continuous linear state feedback. An alternative approach, 
which apparently avoids an appeal to the theory of nonharmonic Fourier series, 
has been offered by Clarke and Williamson [2]. Perturbation theorems for 
spectral operators have been employed by Sun [18] to show that condition (8) 
is both necessary and sufficient for a class of systems that would include the 
one studied in this paper. The main contributions of this paper include a 
more explicit representation of the transformation carrying (2) into (3) than 
that given in [14] as well as a formula for the expansion coefficients of 
in terms of the desired set of eigenvalues 


1. A Discrete Finite Dimensional System 

For the sake of motivation, we will review briefly the reduction of the 
finite dimensional discrete system, where A e 1 ^^^ and b el?, 

\ (l-l) 

to its causal canonical scalar equation 

+ ^l^k-l + •“+ ^nyk-n = '^k* 

where det(Xl -A) = X^+a,X^^ + ***+a. 

1 n 
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For r = l,***,n, we have 

’'k+r “ ^\+r <1*3) 

In particular, if Xq = 0, 

X = [a^ • • • Ab b] 

Well known is the fact that the pair (A,b) is controllable just in case the 
nxn matrix C is invertible. 

With find 


u. 


= Cu. 


(1.4) 


n~l 

^k+n ^l^k+n-1 ^n^k ^p.^m^k+n-m ^n^k 

m=0 


n-1 


n-m 




m=0 


£=1 


n n-£ 

= p(A)x + I I a A 

^ 0 1 A ® 

£=1 m =0 


n-m-£ 


bu 


k+£’ 


(1.5) 


where p(A) = A^ + a^^A^ ^ + ••• + a^I = 0 by the Cayley-Hamilton theorem and 
where we have interchanged the order of summation in obtaining the last line. 
The remaining double sum may be rewritten as 


I 


£=i 


n-£ 


I a A 


n-m-£ 


m=0 


m 


bu 


k-hZ 


CM 


\+l 

^k+n 




where 
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Equation (1«5) now reads 


u 


k+1 




+ 


+ a X, = CM 
n k 


lHc+uJ 


( 1 . 6 ) 


(1.7) 


It is important to note that the coefficients appearing on the left hand 
side of (1.7) are scalars. Assuming the pair (A,b) to be controllable, or 
equivalently, the matrix C to be invertible, we define the variable n by 


■\, 1 - 

It is also important that depends only on k+r. This may be 

demonstrated as follows. Let e^. denote the r*"^ standard unit vector in iP; 
thus 




k+n 


( 1 . 8 ) 


■k,r 


T -1 -1 
= e M C X, . 
r k+n 


According to (1.8), 


m 


k+P,l 


= m“^c 


-1 

^k+p+n 


L\+p,nJ 
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Suppose p+q = r. We wish to show that 


\,r \+p,q 


i.e, that 


T -l.-l T -1--1 

eM C X, .=eM C X, . , 

r k+n q k+p+n 


Since 

^k+p+n ^ \+n ^ ^'^k+A+n ’ 

we must have: 

(i) ej = ejM"^C“UPCM = e^A? 

(ii) e^M“^c“^ 2 bu, . . = 0. 

q uti k+n-£ 

That (i) is true follows from the well-known fact (see [10] or [16], e.g.) 

A 

that A is the companion matrix of A; 
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We now define the sequence {yj^} C 


\+i ■ \,i 


Equation (1.8) now reads 


*^k+l 


'k+n 




k+n' 


Multiplying (1.7) on the left by we find that {yj^} satisfies (1.2). 

Let us now see how linear feedback in (1.1) manifests itself in (1.2). 

If 


“k ■ ' "k-l" 


f = row vector. 


then, by 1 #8, 


u, = f CM 
k 


1 

3 

1 


■^k-n' 

• 

T 

• 

• 

• 

= g 

• 

• 

— 1 
1 


1 

( 


(1.9) 


where > • • • >8j^) • Thus, the closed-loop system 


= (A + bf , 


goes over to the closed-loop system 


y, + (a,-g,)y, .+••• + (a -g )y, = 0. 

1 k-1 n "^k-n 
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Let (not necessarily distinct) eigenvalues of 

A. We are assuming that the entries in A are all real; hence 
is symmetric with respect to the real axis in the complex plane. Given any 
likewise symmetric set of complex numbers easy 

matter to construct a feedback gain f e I? for which the eigenvalues of 
A + bf^ are precisely ^^l»****^n^* Indeed, if 

1 n 

+ c, + ••• + c H n (X-v, ), 

1 " i.i I ’ 

then Cj^ e ^ for 1=1, •••,n and hence, referring to (1.10), we must have 

!=!,•••, n 
T 

where a = (a^,»»*,a^), etc. 

Referring to (1.9), we arrive at the Bass-Gura formula [10] for f: 

T T _1 _1 

f^ = (a-c)^M . (1.11) 

We will show that the same formula holds true under appropriate assumptions 
for the class of infinite dimensional systems mentioned above and to which we 



now turn 
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2* Nonharmonic Fourier Series and Controllability 

In this section, we discuss the nature of the transformation which 
corresponds to the matrix C of the previous section. Before doing so, it is 
first necessary to recall some definitions which concern the scalar input 
system 

x(t) =Ja^(t) + bu(t). (2.1) 

Let generate a strongly-continuous group {5^(t) : t e H} on the 
Hilbert space and denote the action of a conjugate linear functional x 
on Tp ei^by <\p,x>. As with the inner product (•>•) on the bracket 
<•,•> is understood to be conjugate linear in the first argument and linear 
in the second. It is also understood that <*,*> has the property that 
<i(^,x> = (^fX) if X z Let us define the weak solution x(*) to (2.1) 
satisfying x(0) = xq z by requiring 

<'p,x(t)> = <^t)%,x > + / <^t-s)*ip,b>u(s)ds, (2.2) 

0 

for all ill Here, {^t)* ; t e ^ denotes the adjoint group 

generated by with domain snd 6 is understood to represent a 

conjugate linear functional on 

Definition 2.1. [ 4 ] : 

The control distribution element b is said to be admissible with 
respect to if b can be identified with a conjugate linear functional, 
also denoted by fa, for which 

(i) 9(h) 



(ii) for any initial state x(0) z and any locally square 

integrable input u, the weak solution x(t) of (2.1) lies 
in for each t > 0. 

We note that if in (2,1), then the weak solution coincides with 

the mild solution 

t 

x(t) =^(t)x- + / ^(t-s)bu(s)ds. (2.3) 

0 

Since x(t) e for each t > 0 [5], 6 in this case is admissible. 

To make (2.2) somewhat more constructive, we introduce the sequence 
of eigenfunctions of with corresponding eigenvalues will 

assume that each is simple and that forms a Riesz basis for 

Definition 2.2 ([7], [11]): 

A sequence elements of a Hilbert space is said to form a Riesz 

basis for if the exist constants d,D with 0 < d < D < «> such that 

every x e may be expanded uniquely in a series 

k 

with 

d^ I |xj^|2 < 0xl|2 < I |x^|2. (2.4) 

k k 

Equivalently ([7], e.g.), a Riesz basis may be thought of as the image of 
an orthonormal basis under a bounded and boundedly Invertible 


transformation : 
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With every Riesz basis there is associated a unique "dual" Riesz basis 

defined by 

\ 

with the biorthonormality property 
coefficients of x with respect to 

respectively. 

Let us denote by {'1'^} the sequence of eigenfunctions of l.e., 

Suitably normalized, is easily seen to be dual to 

Now let 'J' “ (2.2) and set \(^) “ <’l'j^»x(t)> and bj^ = Then 

oxt t ox(t-s) 

Xj^(t) = e Xj^(O) \ ® u(s)ds, (2.5) 

and hence, weak solutions may be represented as 

X(t) = I Xj^(t)<<)j^, 

k 

with the given by (2.5). If sup|Re < ", admissibility reduces, 

via (2.4), to whether or not, for each t > 0, the sequence 

r \(t-s) 2 

{b, J e u(s)ds] e £ , 

^ 0 

for every u e L^(0,t). 

In the course of reducing (2.1) to canonical form, we will carry out 

manipulations similar to those employed in the previous section and which 


s/ 


The expansion 

^***k^’ are (tPj^,x), (x,<|)j^) 


involve the expression 
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T 

/ ^(T~s)6u(s)ds, 
0 


The whole point of the above discussion is that the latter expression may, 
even if b be Interpreted rigorously if the eigenfunctions of 
form a Riesz basis for and b is admissible with respect to J^: 


T t (A(t-s) 

/^t-s)6u(s)ds = I (b, / e u(s)ds)(}>, 

0 k 0 


In particular, if xq = 0, we may write 


T 

x(T) = / ^(T-s)bu(s)ds H^(T)u. (2.6) 

0 


This is the analogue of (1,4). Conversely, if x the problem of 

p 

constructing a control u e L (0,T) which takes the origin to x at time 
T is characterized by a moment problem: 


^(T)u = X, 


or, with X = I 
k 


T (^(T-s) 

b^ / e u(s)ds = Xj^. (2.7) 


Definition 2.3 ([15]); 

The pair (Q^fb) is said to be approximately, exactly controllable in 
time T if ^(T) : L^(0,T) > is densely, boundedly invertible 
respectively. 
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Referring to (2.7), it is obvious that approximate controllability 
requires b. 0 for all k. But controllability also depends upon 

properties of the sequence of exponentials {e }. This paper is concerned 
with a very special class of _Q^'s, namely those whose spectrum coincides with 
the zero set of a function having the form 


/ \ wr . 

p(o)) = e + a^e 


(oT T a)(T-6) 

• + ••• + a + / a(0)e '• * 

“ & 


d6. 


( 2 . 8 ) 


It is assumed that a^^ e 1^ 1=1, ‘“jm with a^ ^ 0, that a(») eL (0,T) 
is real valued, and that 


Let us write 


0 < 6, < ••• < 0 = T. 

1 m 


p(“) = P„(“) + / a(0)e“^^ d0, 

° 0 


and let denote the zero set of the exponential polynomial p^. It Is 

easy to see that the set {o^} as well as the set must lie in a 

vertical strip of finite width in the complex plane* It has been shown [1] 
using the argument principle that an infinite number of do exist and 

that the number of horizontal strip of fixed height is bounded. 

The prototype in this situation is = 6*^-1 which has for its set 

of zeros “ {kirl}. It can also be shown [17], [19], using the above 

properties of {o^ together with the Fourier transform in the complex 
domain, that the sequence of exponentials {e } forms a Rlesz basis for 
L^(0,T) if, in addition, inf |p^(o, )| > 0. The latter condition ensures the 

k O tC 

bounded invertibility of the map $ discussed after Definition 2.2. One can 

also demonstrate using Rouche's theorem that the sequence of zeros 

2 

of p may be Indexed in such a way that {o)^ - e Z and that 
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inf |p^(“i^)| >0 is sufficient to yield the Riesz basis property for the 




sequence {e }. We omit the somewhat detailed proof of these statements 
[19] and return to their relationship with the notion of controllability. 
Let us denote the sequence biorthonormal to {e } by {q.(*)}: 


/ q,(s)e ds = 


0 


The construction of the sequence {qj^(*)} goes back to the work of Paley and 
Wiener [12] and is carried out as follows. One forms the function 


q (w) = 


p((0) 




and notes that The Paley-Wlener theorem ([12], [13]) then 

A 

asserts that is the Laplace transform of a square integrable function 

^^(•) whose support suppCq^^) S [0,T] : 


q^Cw) = e“® qj^(s)ds. 


Thus, q^ may be represented via the inverse Laplace transform 


= -2¥T q£<“)d“- 


(2.9) 


m • 

The sequence constructed is biorthonormal to {e } since 


Jq q£(s)e ds = 



17 


However, in the next section, we shall derive an explicit formula for each 
without any mention of the Paley-Wlener theory. In- effect, the inverse 
transform (2.9) can be computed explicitly in closed form. 

Assuming bj^ ^ 0 for any k, we may now Invert relationship (2.7) with 
the aid of the sequence (qj^(*)}. Indeed, the control 

»(t) - 

k k 

satisfies (2.7). Here, T is understood to be the quantity appearing in the 
definition (2.8) of p. Thus ^=^(T) has the properties 

q^(T-) - (2.10) 

: <f>k (2.11) 

which will be used throughout the remainder of the paper. 

We conclude this section with a summary of the various assumptions made 
and a brief discussion of their significance. 

Assumption 2.4 ; 

(1) generates a strongly-continuous group t^t) ; t e !R} on a 
Hilbert space 

(ii) The sequence of eigenfunctions forms a Riesz basis 

forj/. 

(iii) The eigenvalues {oj^} of .0^ are the roots of p(tiJ) given by (2.8) 
and p'(Wj^) ^ 0 for any k. 

(iv) The control distribution element 6 is admissible with respect to 

ic 

with coefficients bj^ = ^ 0 for any k. ^ = <^^^). 
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Since finite linear combinations of the are dense in the pair 

under the above assumption, is approximately controllable in time T . 

(c.f. [15]). If inf |bj^| > 0 and inf > 0, exact controllability 

k k 

in time T may be established as well [15]. 

Part (ill) of the above assumption essentially limits the scope of our 
results to the class of linear hyperbolic systems described in the 
introduction. For such systems, p'(a;^) * 0 for all k does not hold true in 
general. We make this hypothesis here only to simplify our presentation - 
notwithstanding a good dose of tedium, it can be shown by introducing 

i 

generalized exponentials (l.e., t'^e ) in the case of multiple roots 
that the canonical form presented in this paper remains valid. 

The case of boundary control (see equation (7)), where b e does 
correspond to an admissible distribution b ([15], Theorem 3.1) 

\ " " <|;J(l)V(l)b, 

and, if no component of b vanishes, approximate conntrollability may be 
established [19]. We do not wish to emphasize the abstract character of the 
above assumptions. However, the development of a canonical form is carried 
out most transparently in this setting. 


3. Reduction to Canonical Form 

In this section, we will derive the transformation which carries the 
system (2) to its canonical form (3), and in order to pursue as closely as 
possible the analogy between the pair (A,b) of the discrete system above 
with the pair (o/,b) at hand, we will proceed somewhat formally. Assumption 
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2.4 will be understood to be in force throughout. We begin with an analogue 
of the Cayley-Hamilton theory: 


Lemma 3.1: Under (i), (ii), and (iii) of Assumption 2.4 , 


^(t) + a,5f(t-0,) + ••• + aMt-T) + / a(0)^(t-0)d0 = 0. 
11 m Q 


Proof . Given an arbitrary initial state x = I the solution to 

k 

x(t) =Q^(t) is given formally by x(t) =^(t)x and concretely by 


w, t 

x(t) = I e 
k 


Thus 


Wit) + a^(t-0p + ••• + a^t-T) + / a(0Mt-0)d0]x 

<\(t-T) 

= I P(<\)e = 0. 

k 


This completes the proof of the lemma. 


Formally, solutions to x(t) =J^(t) + bu(t) satisfy 


x(t+r) =if(r)x(t) + /^(r-s)bu(s)ds. 

0 


(3.1) 


Thus , 
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T 

x(t+T) + a,x(t+T-e.) + ••• + a x(t) + / a(0)x(t+T-e)d0 
1 1 m A . 


T 

= l^(t) + a^t-ep + ••• + a^(t-T) + / a(0)i^t-0)d0]x(T) 


+ / ^(T—s )6u(t+s )ds + a./ ^(T— 0. — s)bu(t+s)ds 

0 0 


+ ••• + a , / -^(T-0 ,-s)'bu(t+s)ds 

m— 1 * m—i 


T T-0 

+ I a(0) / ^(T-0-s)6u(t+s)ds d0. 

0 0 


(3.2) 


The first term Is zero by Lemma 3.1. Changing the inner variable of 
integration from s to s+0 and then interchanging the order of integration, 
the last term may be written as 


T s 

/ ^Cr-s)b f a(0)u(t+s-0)d0 ds. 
0 0 


Likewise, 


/ ^(T-0.-s)bu(t+s)ds = /if(T-0 -s)bxfo T-e i(s)u(t+s)ds 

0 0 ^ ’ i^ 


T 

= / ^(T-s)bXjQ^,p_0 j (s-0^)u(t+s-0^)ds 


T 

= f ^(T-s)bx^Q (s)u(t+s-0^)ds. 


where characteristic function of the interval [a,b]. The 

right-hand side of (3.2) may now be written as 
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/ ^(T-s)b(^u(t+0)(s)ds H^u(t+*), 
0 


where So is as before and JC* ^ (0,T) > L"^(0,T) is the transformation 
defined by 


C^y)(s) = y(s) + ajX [0 (s)y(s- 0 j) + ••• 


i^rfi *'■ / a(0)y(s“0)de. 


(3.3) 


Equation (3.2) now reads 


x(t+T) + a,x(t+T-0,) + ••• + a x(t) + / a(O)x(t+T-6)d0 =^^u(t+*). (3.4) 

i i m Q 


Equations (3.3) and (3.4) should be compared with (1*6) and (1.7) 
respectively. To proceed further, the nature of the transformation ^ must 
be analyzed. 


Lemma 3.2: 

( i ) ^ is bounded and boundedly invertible . 

(ii) If g e L^(0,T) with supp(g)c [T-0,T] for some 0 e (0,T), then 
suppj^ ^g)C [T-0,T] as well. 

Proof: That ^ is bounded is clear from (3.3). Referring to the 

latter, let us write 

s 

with C^y)(s) = - / a(0)y(s-0)d0. A standard theorem from integral 

^ 0 
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equations [8] yields the fact that exists and is bounded with 

00 

^ k=l 

It is easily established that suppC/^g) and hence ^g) £ [T-0,T] 

if supp(g) £ [T-0,T]. Moreover,^ is nilpotent since 
suppOf'g)® [T-0+0j,T] if supp(g) £ [T-0,T]. Thus is 

nilpotent and 


= c^-4)"^0+I 

k=l ^ 

for some n < «, which shows that is bounded and satisfies (ii). This 

completes the proof • 


We remark that part (ii) of Lemma 3.2 is the analogue of the fact that 

span{e ,***,e } is invariant under the matrix discussed in Section 

n— p n 

01 • 

1. We next Fourier analyze the action of ^ on the exponentials e : 


Lemma 3.3: 


/ )(s)ds = p(o))/( 




V' 


Proof : 


In forming the expression on the left, we find terms of the form 


T 

/ e 
0 


^[0.,T] 


ds , 


and the term 


T o^(T-s) 
/ e 
0 


s w (s-0) 

/ a(0)e^ 

0 


d0 ds . 



Standard manipulations reduce the. former to 


and the latter to 


o)(T-0^) 

e 


e 


( 0 - 0 ), 

k 


> 


I a(9) ® 
0 


a)(T-0) Wj^(T-0) 


- e 
0) - u. 


d0. 


It follows that 


J e X#e )(s)ds = [e - e + aJ e 


- e 


+#»«+a,le -e 1 

m-l'^ ^ 


T aj(T-e) 01 (T-0) 

+ / a(0)(e - e )d0]/(o)-a) ) 

0 


= (p(“)-p(«^))/(“-\) = p(a))/(aj-(Uj^), 


since o)^ is a zero of p. This completes the proof. 

2 

Recall that the basis {e } CZ L (0,T) has a unique biorthonormal 

2 

sequence which we denote by (0,T). The key property which 

possesses is the following: 


Lemma 3.4: 


^ : e ^ p'(o)^)q^(T-). 
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Proof : According to Leimna 3*3, 


T u) ox. • 

/ e ^ Cde )(s)ds = 


Since the sequence 
conclude that 


(qj(T-)} 


Is biorthonormal to 


0) (T-«) 
{e J } 


and unique, we 




\ 


)(s) = p'(‘>ij^)qj^(T-s), 


which completes the proof. 

Remark ; 

The preceding lemma allows us to construct the biorthonormal sequence 
explicitly: 

1 

q^Cs) )(T-s). 

Ic 


Under Assumption 2*4, (iv), the map 5^ possesses a densely defined 
inverse (see (2.11)). If x(») is a solution of x + bu with 

sufficiently regular initial state xq and input u(*), we may define 


n(t,*) = I<«r"^"^x(t+T))(*) e L^(0,T). 


(3.5) 


Let us assume that n Is defined for every t and s. Then 


Lemma 3.5: n(t,s) is a function of s+t. 


Proof : Suppose p+s = r < T. We wish to show that n(t,r) = n(t-hp,s). 


According to (3.5) , 
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n(t,r) ^^■^’■^xCt+DjCr). 


If we let x(t+T) = 'I Xj^(t+T)<j)j^ and use (2.11) and Lemma 3.4, we have 


-1 -1 

n(t,r) = I p'(w, ) b. X. (t+T)e 


Likewise, 


k' ‘'k "k' 


n(t+p,s) = I P'(\) ^bj^^Xj^(t+p+T)e 
k 




Since 


01 p p 01 (p-0) 

Xj^(t4p+T) = e Xj^(t+T) + bj^ / e u(t+9+T)d0, 


we have 


:.r) = I (p'(oi ) M e u 


n(t+p,s) - ri(t,r) = 2. Ip'C^^) 

k 0 


01 s 

(t+0+T)d0)e . 


If we change the variable of Integration via p-0 = T-x and invoke Lemma 3.4, 
the right-hand side becomes 


-1 T (T-t) _ 

= ^ I / e u(t+p+x)dT q (T-»)}(s). 

k T-p 


But the function on which operates is nothing but 

expanded in terms of the sequence {q^^CT-* )} . Hence 


^[T-p,T] 


( •)u(t+p+*) 


= C^"^g)(s). 


Clearly, supp(g) <z [T-p,T], and since 0 < p+s = r < T implies s e [0,T-p], 

\ 
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we conclude from Lemma 3,2 that X^^g)(s) = 0. Thus, n(t4p,s) = n(t,r)* 

This completes the proof. 

Lemma 3.5 may be shown to hold true in the case where n(t,*) is merely 
a square-integrable function for each t by approximating everything with 
continuous functions and then passing to the limit. We may now define the 
scalar variable 

y(t+0 = n(t,‘) = U"V’^x(t+T))(*). (3.6) 

We conclude this section with 

Theorem 3.6: Under Assumption (2.4), the variable y defined by 

(3.6) satisfies the functional equation 

T 

y(t) + a y(t-0.) + ••• + a y(t-T) + / a(0)y(t-0)d0 = u(t). (3.7) 

11 m Q 

Proof; Apply on the left to both sides of (3.4). 


4. Spectral Synthesis 

Let us first examine how continuous state feedback 

u(t) = ((J,x(t) 

in the original system (2.1) manifests itself. in the canonical system (3.7). 
Let have the expansion 



where 


{tj} denotes the sequence of eigenfunctions of 
y(t+») have the expansions 


Let x(t) and 


X(t) = I Xj^(t)(|)j^ 

k 

(4.1) 

y(t+*) = I yi,(t)e ^ . 
k 


Referring to (3.6), we have 


Xj^(t+T) = \p'(\)yk(t). 


Thus, 

u(t) = (i5,x(t)) = I fk\p'(\)yk(*^"’^)* 

k 

Again denoting the sequence biorthonormal to {e } by {qj(*))> the right- 
hand side of the last equation may be written as 


/ (I f.b.p'(“.)q.(s)]( I y. (t-T)e'*^ )ds. 
0 j j J i i j k 


By (4.1), 


w s 

I y. (t-T)e ^ = y(t-T+s), 

k 


If we change the variable of integration via s = T - 0 and define 
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we obtain 


T 


u(t) = / 

0 


g(0)y(t-e)de. 


The resulting closed-loop canonical system Is simply 

T 

y(t) + a,y(t-6,) + + a y(t-T) + / (a(0)-g(6))y(t-6)d6 = 0. (4.3) 

11m n 


Thus, the eigenvalues of the closed-loop system 


x(t) “ (j^ + 6 <8> (( )x(t) 

are roots of 

T 

p (m) = p ( 10 ) + / (a(0) - g(0))e“^'^"®^d0. (4.4) 

® ° 0 

A-- 

As we have already mentioned, in Section 2, one can show using Rouche's 

theorem that the sequence {\) of roots of Pg differs from the sequence 

2 

{a^} of roots of p^ in the sense that e Z • This suggests that 

any sequence with this property can be synthesized with a certain 

feedback gain ejk 


Lemma 4.1: is a sequence of distinct complex numbers for 

2 ' 2 
which - 0 ^} e £ , then there exists a unique c(*) eL (0,T) for which 

the zero set of 

T 

P («o) = p (m) + / c(0)e“^'^"®^d0 

° 0 

coincides with { ^ • 
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Proof? c is characterized by the moment problem 

T . V (T-9) 

p^(Vj^) = 0 ==> / c(6)e d0 = 


Using the properties of the exponential polynomial p^ discussed in Section 2 

2 

together with the assumption that e A j it is readily established 

2 

that {p (v,)} e £ • Moreover, it has been shown [17] that the sequence 

O K. 

V • 

k 9 

{e } forms a basis for L^(0,T). Denoting the sequence biorthonormal to 

\* 

{e } by {hj^(*)}, the unique solution to the above moment problem is given 
by 


c(-) = -I Po(\)\(T-), 

which completes the proof. 

We now come to our main theorem. 


Theorem 4.2: 
for which 


Let { v^} be any sequence of distinct complex numbers 


k - 


e £ 


and let Assumption 2.4 hold true. Then there exists a unique ^ for which 
the spectrum of b <8 coincides with { \ } • 


Proof . By Lemma 4.1, we may express the desired closed-loop 

characteristic function uniquely as 


T 

p^(m) = p (w) + / 

° 0 


c(0)e 


o)(T-0) 


d0. 


Referring to (4.2) and (4.4), the desired feedback gain element 
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i = I 




must be such that 


P„(“) = P„(“) + / (a(0) - g(6)}e“^'^ = p (w), 

0 ^ 


or g = a - c. Let a and c have the expansions 




:(•) = I Cj^qk(T-) 


Referring to (4,2), we must have 


\ - '=k> 


or 


*k ■ '“k - <=k>p'<“k>'‘'’k‘ 


(4.5) 


This should be compared with (1.11). It remains only to show that our 

2 

hypotheses ensure that ^ ^ expansion coefficients aj^ and Ck 


T \ (T-6) 

\ = / a(0)e d0 = p(<^) - Pq(\) = 


T to, (T-0) 

Cy. = I ^ c(0)e d0 = P^(\) - Pq(\)' 


are given by 
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Thus, by the intermediate value theorem, 


- c, 


= -Pc^V = -Pc<V - 


for some on the line segment connecting and 


sup|p'(\)| 


< «», we conclude that 


Since 


{■ 


\ \ I 

P'(\V 


e r 


just in case 


\ 


} e I 


2 


This completes the proof. 


For linear hyperbolic systems of the type described in the introduction, 
is the real operator 

-Q^= A(x) ^ + A(x), 

with domain 

9(0> = = 0 = <!.■*■( 1) - D^<|)”(1)}, 

where is the Sobolev space of square integrable <f^’^-valued functions 

possessing a square integrable derivative. As such the eigenfunctions of 
as well as those of will have the property that Likewise, if 

{V|^} is a sequence of complex numbers with ^ k* biorthogonal 

set of functions (h^(*)} associated with the exponentials {e } will have 
the property that h^ = ^-k* Thus, for any such sequence {^} satisfying 
the hypothesis of Lemma 4.1, the function c(*) will be real-valued with 
expansion coefficients satisfying c^ = hence the feedback gain 


element 
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k 


= I 

k 


a, - e, 


> 


constructed above will have each of its 2n component functions real valued. 

We conclude this section with a brief discussion of how to use the above 

theory to construct a feedback gain which shifts a finite number of 

eigenvalues. Suppose we wish to shift o) ,03 , to ^ij***,^ • The 

i z n 1 n 

desired closed-loop characteristic function is thus 


p^(u) 


g(^) 

3(w) 




n n 

where a(o)) = IT (w-v ) and B(o)) = Tl (co-o) ) . Thus 
k-1 ^ k=l ^ 


Pc(0)) 


p(to) + 


g((u) - g(g)) 
S(io) 


p(u) . 


Since the degree of a - 6 Is strictly less than n, (a-0)/S admits a 

t 

partial fractions decomposition 


g(m) - 6(m) _ ^ 

B(») " til « - \ • 


(4.6) 


According to Lemma 3.3, 


i U3 • 

*d0. 


0) - OJ, 


Thus 

T n 10 • ox 

p (“) = p(to) + / I r.e )(0)e^ ~ ^d6 . 

0 k=l 


By Lemma 3.4, this is 
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T T1 

,((o) = p(a>) + / I r.p'(03. )q. (T-0)e“^^"®^d0 
0 k=l ^ 


= p(0)) + / c(0)e“^^“®^d0. 
0 


Thus, the desired expansion coefficients of c are given simply by 
c^ = with the residues r|^ being determined by (4*6) • 


5. Concluding Remarks 

In this paper, we have demonstrated precisely to what extent the 
eigenvalues associated with a controllable pair 6) of a certain type may 

be modified via continous linear state feedback. Our results parallel those 
of Russell [14] for the case m = 1 as explained in the introduction. Our 
main contributions here include a more detailed analysis of the canonical 
transformation and a direct method for computing feedback gains in the case of 
shifting a finite number of eigenvalues. All of this requires, of course, a 
rather detailed knowledge of the spectral structure of In practice, one 

can only obtain this information approximately by numerical computation. In 
this regard, our results might prove useful in providing exact solutions in 
special cases which could then be used to determine the accuracy of numerical 
computations in more general cases. It should be pointed out that it is 
unrealistic to assume the availability of state feeback in systems of the type 
discussed in this paper. Nevertheless, we feel that the results herein 
obtained should play a role in the development of an observer or asymptotic 
state estimator theory for such systems. 
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